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Abstract 

We develop an analytic method and approximations to compute the 
polarization induced in the cosmic microwave background radiation on 
a wide range of angular scales by anisotropic Thomson scattering in 
presence of adiabatic scalar (energy-density) linear fluctuations. The 
formalism is an extension to the polarized case of the analytic approach 
recently developed by Hu and Sugiyama to evaluate the (unpolarized) 
temperature correlation function. The analytic approach helps to high- 
light the dependence of potentially measurable polarization properties of 
the cosmic microwave background radiation upon various parameters of 
the cosmological model. We show, for instance, that the ratio between 
the multipoles of the temperature and polarization correlation functions 
depends very sensitively upon the value of Qo, the matter density in 
units of the critical, in an open universe. 

PACS numbers: 98.80.-k, 98.70.Vc, 98.80.Es 
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I. INTRODUCTION 



The cosmic microwave background radiation (CMB) is a seemingly unlimited source of 
information about the history and evolution of the Universe. Its existence is one of the most 
powerful arguments for the standard "hot" big-bang cosmology. Its high degree of isotropy 
is an indication of the Universe large scale homogeneity. Its precise black-body spectrum 
places stringent bounds on alternative cosmologies. Its large scale anisotropies uncover tiny 
fluctuations in the gravitational potential, most likely the same that once upon a time led 
to the formation of galaxies and other large scale structure in the Universe. 

The polarization properties of the cosmic microwave background radiation constitute 
yet another set of observables whose eventual measurement is still to bear fruit. A positive 
measurement of a degree of polarization would provide much further insight into the Universe 
history and evolution. The current bound on the degree of linear polarization of the CMB 
on large angular scales, P < 6 x 1CT 5 JT|, already serves to set limits and constraints upon 
alternative cosmological scenarios. 

The fact that anisotropic Thomson scattering of photons and electrons around the time 
of decoupling induces a degree of linear polarization in the cosmic microwave background 
was pointed out shortly after its discovery in 1965. Several estimates were made of the 
degree of polarization expected. The task is to evaluate the Stokes parameters of the pho- 
ton distribution function, which satisfies a Boltzmann equation with a Thomson scattering 
collisional term |J. The first computations were performed in homogeneous but anisotrop- 
ically expanding universes Energy- density fluctuations P-|15|1 and long wavelength 

gravitational waves [|3~6| , [T7l , n~2| , p~3| , p~8f| were also considered as the source of the anisotropy that 
leads to polarization. Typically, at least within standard recombination histories, the degree 
of linear polarization turns out to be almost two orders of magnitude smaller than the large 
scale temperature anisotropy. This is still well beyond current detection capabilities. Never- 
theless, with large scale anisotropies in the CMB now positively measured by COBE-DMR 
|T9|| and other experiments, it is worthwhile to further analyse the predictions for the CMB 
polarization properties within alternative realistic cosmological models, and its dependence 
upon various parameters. 

Very detailed numerical computations of the CMB temperature fluctuations and polar- 
ization correlation functions on all angular scales were performed by Bond and Efstathiou 
lT| for standard Qq = 1 cold dark matter dominated universes, with scale invariant adia- 



batic or isocurvature scalar fluctuations. It is one of the purposes of the present paper to 
reproduce the most significant features of the multipole expansion of the polarization corre- 
lation function on all angular scales in this cosmological model from an analytic approach, 
to better highlight its dependence upon various parameters. To achieve that goal, in Section 
II we develop a method and perform approximations that very closely parallel the recent 
analytic approach introduced by Hu and Sugiyama to evaluate the temperature anisotropies 
on all angular scales for the unpolarized case |2(| . In Section III the multipole expansion of 
the polarization correlation function is performed, and an expression for the total degree of 
linear polarization is derived. In Section IV we specialize to a specific cosmological model, 
namely a spatially- flat, Qq — 1 cold dark matter cosmology with scale invariant scalar fluc- 
tuations normalized to the large angle anisotropy measured by COBE, under the assumption 
that the baryonic energy density is much smaller than the radiation energy density at the 
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time of decoupling. In Section V we find the functional dependence of the ratio between the 
multipoles of the temperature fluctuation and of the polarization correlation functions upon 
Qq, the mass-density in units of the critical, both in open universes as well as in spatially- flat 
models with a non- vanishing cosmological constant, with scale-invariant scalar fluctuations 
as the source of the anisotropy. Section VI rounds up the conclusions. 



II. ANALYTIC APPROACH 

A. Tight Coupling approximation 

The Boltzmann equations for the photon distribution function in realistic cold dark 
matter cosmological models with energy-density linear fluctuations have been numerically 
solved, and the predictions for the CMB temperature correlation function on all angular 
scales have been analised in great detail j9|, plljriH . The Boltzmann code was also used to 



study the dependence of the microwave background anisotropies upon various cosmological 
parameters |22|. Notwithstanding the thoroughness of these numerical studies, it is very 



useful to have some analytic approximation to the exact solution, to gain even further 
insight into the nature and origin of the CMB anisotropies. Some analytic methods and 
approximations to compute with reasonable accuracy the CMB anisotropies on all angular 
scales were recently developed p0|j23|j2^| , improving upon the original work of Sachs and 



Wolfe m 



Here we extend to the polarized case the analytic approach that Hu and Sugiyama 
introduced to compute CMB anisotropies J20j . The method is in turn an extension of the 
standard tight coupling analysis p6| , p7| to include realistic time dependent gravitational 



potentials, described in a gauge invariant formalism |]28|- p0|1 . It is based upon an expansion 



of the CMB temperature fluctuation in inverse powers of the differential optical depth. In 
the tight coupling regime, when the effectiveness of Thomson scattering of CMB photons 
with free electrons makes the differential optical depth high, a perturbative expansion to first 
order constitutes a very good approximation to the exact result. While photons and baryons 
are tightly coupled, all higher multipoles of the temperature fluctuation can be evaluated in 
terms of the monopole, which in turn obeys the equation of a forced and damped oscillator 
20fl . We will see now how to include the polarization dependence of Thomson scattering 



into this formalism. 

The CMB gauge-invariant temperature fluctuation for a given direction of observation, 
described in terms of polar angles 9, <fi, is the relative temperature fluctuation around the 
mean evaluated in the shear-free (Newtonian) gauge: A T (9,4>) = AT(9,(f))/T. We compute 
first the temperature fluctuations induced by only one Fourier mode of the scalar fluctuations 
in the gravitational potential, with wave-vector k. Then it is convenient to use a reference 
frame such that z \\ k, since there is axial symmetry around the direction of k. The degree of 
linear polarization, Ap, is defined in terms of the Stokes parameters Q and U @ of the CMB 
radiation. We choose the two orthogonal directions into which the intensity is projected to 
define the Stokes parameters as 9 and <fi of spherical coordinates. The advantage is that, in 
this basis, U = 0, and the only non- vanishing Stokes parameter is Q = I§ — l^j Iq + I^ = Ap 
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Given the axial symmetry around the direction of k \\ z, the multipole expansion of either 
the temperature fluctuation or the polarization in a given direction of observation n can be 
written as 

A(n,£) = £(2/ + l)A^(/i) (2.1) 
i 



with Pi the Legendre polynomials and [i = cos 6 = k ■ n/k |32[ 



In a spatially-flat Robertson- Walker metric with linear density fluctuations described 
by gauge invariant potentials $, and after angular integration of the collisional term of 
the Boltzmann equations 0, made easier by the axial symmetry, the evolution equations 
for the Fourier mode of wavevector k of the gauge-invariant temperature fluctuation and 



polarization read [20.11.33 



A T + ikfi(A T + = -6 - k{A T - A T0 - yVb ~ |P 2 (/i)[A T 2 + A P2 - A P0 ]} 

(2.2) 

A P + ik/jAp = -k{A P + |[1 - P 2 (y)][A T2 + A P2 - Ap ]} 

where a dot means derivative with respect to the conformal time r = / dta a /a, with a(t) 
the scale factor of the spatially-flat Robert son- Walker metric, and a Q = a(t Q ) its value at 
present time. We shall write in this work the present value of the Hubble coefficient as 
Hq = /ilOOkm.s - .Mpc~ . k = x e n e GTCL/ 'clq is the differential optical depth for Thomson 
scattering, with x e the fraction of ionized electrons with number density n e , and <jp the 
Thomson scattering cross section. V& is the velocity of the baryons. $, ^ are the Fourier 
modes of the gravitational potentials [^8 ,[30 . 



The equation of motion for the baryons reads 

V b = -~V b - + ^(3A T1 - V b ) (2.3) 
a K 

where R = 3p&/4p 7 , the ratio between baryonic and radiation densitites. 
Equations (|2.2j) can be formally integrated 



(2.4) 



(A T + *) = J ro dTe ik ^ T ~ To) e-< TQ ^ 

{k(A T0 + * + fiV b + |P 2 (^)[At2 + A P2 - Apo]) - 6 + ^} 



jr dre iMr-7D) Ae -*(70,T)I[ 1 _ P 2 ^)][^ T2 + A P2 - A P0 ] 



where 



I a ( T ) 

k(t ,t) = / x e n e a T ——dT (2.5) 
J T a(r ) 

is the optical depth to photons emitted at conformal time r. The combination ke~ K is called 
the conformal time visibility function. It is the probability that photons last scattered within 
dr of r. For standard recombination this function has a sharp peak at the conformal time 
of decoupling tjj p4]| . Thus, the integral for Ap in eq. ( |2.4j ) is dominated by the value of 
the integrand around decoupling. In other words, for standard recombination histories, with 
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no reionization, the polarization of the CMB we observe today was produced just before 
decoupling. 

In order to approximately solve the integral for Ap in eq. ( |2.4| ), we need to know the 
value of the combination Sp = [A PQ — A T2 — Ap 2 ] around decoupling. We can evaluate the 
first multipoles of the CMB temperature fluctuations and polarization before decoupling in 



the tight coupling approximation [20], i.e. by a perturbative expansion in inverse powers of 
the differential optical depth k, which is high before decoupling. In order to perform this 
perturbative expansion we first rewrite the evolution equations ( |2.2| ) in the following form 



At - A T0 - fiV b - |P 2 (/i)(A T2 + A P2 - A P0 ) = -r c [A T + ikfi{A T + ) + $] 

A P + |(1 - P 2 (/i))[A T2 + A P2 - Ap ] = -T C {A P + ikfxAp) (2.6) 

3A T1 -V b = r c R[a~ l £(aV b ) + ik% 

and expand all quantities of interest in powers of re = kT 1 , the conformal time between 
Compton (or Thomson) scatterings, assumed large for the tight coupling to be a sensible 
approximation. 

In the strict tight coupling limit, that is to order zero in tq = ft" 1 , the solutions to these 
equations read 



A Tl = \V b ; A TZ = if l>2 
A P = . 



(2.7) 



The interpretation of these formulae is very simple. In the lowest order approximation the 
photons and baryons are so strongly coupled that the photon distribution is isotropic in the 
baryon's rest frame. The photon distribution being isotropic, Thomson scattering does not 
polarize the CMB. 

Now we expand eqs. (|2.6 ) to first order in tq = k~ l , and get 



Ap 2 = -|A P0 = iA T2 ; A T2 = --^ikr c A T1 

A T1 = f(A ro + <i>) (2.8) 

A Tl = A Pl = if I > 3 ; A P1 = . 

These equations also have a simple interpretation. The polarization of the CMB is pro- 
portional to the quadrupole of the photon distribution function (a dipole does not induce 
polarization). The quadrupole in the temperature fluctuation, in its turn, is produced 
by the "free streaming" of the dipole between collisions. We see this from the relation 
A T2 kr c A T1 . The tight coupling approximation is actually valid when kr c « 1, i.e. for 
wavelengths much larger than the photon mean free path. The dipole in the temperature 
fluctuation can be derived from the gravitational potential and the monopole, through the 
relation ikAn = — {Apo + The monopole itself obeys, in the tight coupling limit, the 
equation of a forced and damped oscillator |20| 
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x i? • k 2 ■■ R ■ k 2 ^ . , 

The solution to this equation, for a given cosmological model, explains the properties of the 
CMB anisotropies ||20|| , and will also determine its polarization properties, as we show in the 



next subsection. We postpone to Section IV the solution of eq. ( |2.9|) for a specific, realistic 
cosmological model, and the analysis of its most significant features, such as Doppler peaks, 
etc. 



B. Polarization 

Here we evaluate the integral for the polarization Ap in eq. fl2.4| ), which can be rewritten 

as: 

where we have defined 



Ap = -(l-/i 2 )/ dre ikll{T ~ To) ke- KM Sp (2.10) 
Jo 



S P = Ap — A T2 — A P2 . (2.11) 

Because the visibility function ke~ K ^ r>,T ' is strongly peaked around the time of decoupling, 
t d , to calculate the polarization today it is only necessary to know Sp near decoupling. 
Right before decoupling, and to first order in the tight coupling expansion, valid for scales 
such that krc « 1, Sp can be approximated, using eqs. (|2.8|), as 



Spn^ihrcATi ■ (2.12) 

But since tq = k,^ 1 grows very fast during recombination, we need to know the time de- 
pendence of Sp around decoupling with better approximation than this. From (|2.2|) we find 
that Sp satisfies the following equation 

S P + ^-kS P = ik[-A T1 + -(A T3 + Ap 3 - A P1 )] (2.13) 
10 5 5 

During recombination we can approximate the right hand side of this equation by its tight 
coupling expansion to first order, and thus keep just A T1 , as given by eq. ( |2.8| ). In this case 
the approximate solution for Sp is 

S p (t) = -ik f r dr'A T1 e-ro^') (2.14) 
5 Jo 

We now replace this in the integrand of eq. ( |2.10|) for Ap. The integral is dominated by 
the contribution around decoupling, since the visibility function is strongly peaked around 
Tp>, the conformal time of decoupling. The conformal time visibility function is well approx- 
imated by a gaussian, of width Arp> p5|,|34]]. This means that photons were able to travel a 



distance of order Atd between their last two scatterings. This is the time the quadrupole 
had to grow, and thus the final polarization should be proportional to /cAt^Api. To see 



6 



that this is indeed the case, we perform the integrals leading to Sp(r) around decoupling 
and to Ap under the following approximations, analogous to those in refs. [16|,13l|. We first 



approximate «(to,t) ~ — , which is justified by the gaussian nature of the visibility 
function during recombination. Notice also that k(t,t') = k(t ,t') — k(t ,t). We also ne- 
glect the time variation of An during the decoupling transition. Then we can write, for r 
around decoupling: 

9 o , . POO rlqr , 

S p (t) « -tkA T1 (r D )Ar D e^,r) / _± e -&- ? (2.15) 
5 Ji x 

where we have changed the integration variable from r' to x — k(t ,t)/k(t ,t'). Now, 
neglecting also the time variation of e ik p( T - T °) during recombination, we get 



A P = |(1 - /j 2 )e ik ^ TD - T ^hkA T1 (T D )AT D / °° dne~To K Jj~ -< ™ 



oo dx „ — xk 

(2.16) 

(1 - /i 2 )e ifc ^ TD - r °)0.51iA;A r i(rp)Arp, = (1 - ^ 2 )e ik ^ TD ~ T ^ (5{k) . 



We here defined, for shortness and later reference, the quantity (3{k). Expression Q2.16 ) 



is one of our main analytic results. It gives, for standard recombination histories, the 
polarization induced upon the CMB by one Fourier mode of wavevector k of the linear 
density fluctuations, in terms of the value of the dipole in the total temperature fluctuation 
at the time of decoupling. It is proportional to the width of the last scattering surface, 
Atd, because it is actually the quadrupole in the temperature fluctuation during the last 
few scatterings what induces the polarization, and the quadrupole itself is proportional to 
the dipole times the width of the last scattering surface. 

Expression ( [2.16|) is strictly valid only for scales such that kArp << 1, since we took 



the exponential out of the integral and simply evaluated it at r = Tp. For scales such that 
kArp) » 1, the oscillations in the integrand produce a cancellation. In other words, the 
finite thickness of the last scattering surface damps the final polarization on these scales. 
The degree of polarization is thus largest for modes of wavelength comparable to Arp. 



In the long wavelength limit, the expression ( 2.16 ) reduces to our previous analytic 



estimate of the polarization on large angular scales |13| , based on the method developed by 



Polnarev fl6| , once the dependence of Api with the gravitational potentials is replaced. 



The polarization properties of the CMB are very sensitive to the details of the ionization 
history, and could very well serve to trace it back P,p6| . The proportionality in Atd in our 
solution is a hint of this. In a scenario with an appreciable reionization, the polarization 
would increase, because the quadrupole of the temperature anisotropy in the electron's rest 
frame, which is the source of polarization, would be greatly enhanced. 

A very important conclusion that can be drawn from eq. ( 2.16|) is that one can deter- 



mine the value of the dipole of the temperature fluctuations at recombination measuring 
the present polarization properties of the CMB radiation, at least if there was no reioniza- 
tion after recombination. This is a very interesting perspective, since it could also serve, 
for instance, to test alternative evolutions after recombination. The value of the dipole at 
recombination, Api(rp), can be derived from the monopole and the gravitational potentials 
using eq. (^.8[), and the monopole itself solving eq. (|2.9|), in the tight coupling approxima- 
tion. The oscillatory behaviour of Ap corresponds to the so-called Doppler peaks. In the 
case of adiabatic density fluctuations, the monopole turns out to be proportional to cos</> 
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with <fi — k Jq drc s , where c s = [3(1 + i?)] -1 / 2 is the photon-baryon fluid sound speed, while 
the dipole is proportional to sin <fi |[20fl . Thus, in the case of adiabatic perturbations the peaks 
in the polarization Ap are located at wavevectors such that <P{td) = (m+ with m an 
integer. For models with low ilb where c s ~ const. ~ 4=, the peaks are at = (m+|)7r. In 
the case of isocurvature perturbations the monopole is proportional to sin (ft, and the peaks 
in the polarization are instead at 4>{td) = mir. A test that the relative locations of the 
peaks in At and Ap verify these relations may serve as a test if the recombination process 
was the standard one assumed here or not. 

The relative heights of the different peaks is also a potential source of information about 
cosmological parameters. As can be seen from eq. (|2.4j ), when Ap and \1/ have opposite 
sign a supression in the height of the corresponding peak in A<r may occur. In the case of 
adiabatic fluctuations, this can happen to even peaks. This pattern of suppresion of even 
peaks relative to odd ones is sensitive to the value of VL b h 2 : the suppression grows with VL b h 2 . 
pOfl This pattern of relative suppresion in the heights of even peaks does not occur to the 
polarization. This is so because the gravitational infall represented by which affects the 
anisotropy, does not change the degree of polarization. 

The CMB polarization is proportional to the dipole in the temperature fluctuation at 
recombination, which in turn, being proportional to a time derivative of the monopole, 
is proportional to c s , in the region where A T0 oc cos (ft, with (f> = k J T drc s . Since c s oc 
(1 + R)~ l l 2 and R oc h 2 Clb, then the height of the peaks in the polarization decreases for 
larger Qf, (any other dependence on Q b in eq. ( 2.16|) is not very significant). 

More conclusions will be drawn in Sections IV and V, in the context of more specific 
cosmological models. Now we turn our attention to the behaviour of the temperature and 
polarization fluctuations at smaller scales, where the tight coupling approximation starts to 
break down. 



C. Diffusion damping 



The approximations that lead to eq. ( |2.9| ) break down for scales much smaller than Tc 
(krc » 1), when the coupling between photons and electrons is not so tight. To find 
the qualitative, and approximately quantitative, behaviour at small scales it is necessary to 
expand the temperature, polarization and velocity fluctuations to second order in tc = k . 
For these scales the role of gravity is unimportant, so we solve equation ( |2.2| ) neglecting the 
gravitational potentials |]20|j37| . Assuming solutions of the form 



A T (r) = A T e^ , A P (r) = A P e^ , V b (r) = V b e WT , (2.17) 



substituting this ansatz into the evolution equations ( |2.6| ), and expanding to second order 
in rc, we obtain for u = uq + 27 



where we have defined c s , the photon-baryon sound speed. 



(2.18) 



S 



The most important lesson to derive from this result is that at small scales, those 
such that k >> kp, defined below, the CMB temperature fluctuations and polarization 
are damped by an exponential factor e~^, with 

k 2 , 9 f T dr 1 r 9 16 . . 

^?^/,Ts(iW |S+ B (1 + fi)l ' (2 ' 19) 

where we have taken into account the evolution of R and k with r. This is just Silk damping 
due to photon diffusion [|38|, |37| , f2"0H . 



It is worth to stress at this point that the rigorous derivation of the damping factor 7 
requires to take into account the polarization dependence of Thomson scattering, as we did 
here and was first pointed out by N. Kaiser |39] . If the evolution equations for Ap are solved 
neglecting polarization (ignoring Ap in eqs. (|2.2j )), then a factor 4/5 would appear instead of 
the factor 16/15 in expression ( |2.19|) for 7 |37].|2(|. The polarization dependence is important 



because after each scattering the radiation is partially polarized. The polarization acts as a 
source of anisotropy and the anisotropy as a source of polarization. This coupling between 
anisotropy and polarization makes the fluctuations decay more rapidly than if there were 
no polarization dependence |J5|. The difference is not insignificant. For a CDM dominated 



universe such that R -C 1, before recombination (x e = 1) we get 

k D 2 « 10 7 (1 - Y P /2)-\tt b h)- l (l + z)~ 5 / 2 h- 2 Mpc 2 (2.20) 
which differs by about 20% from the result obtained neglecting that the radiation is polarized 



20. 



III. CORRELATION FUNCTION AND TOTAL POLARIZATION 

In this section we perform the multipole expansion of the CMB polarization correlation 
function, and evaluate the total polarization produced by the scalar fluctuations. In the 
previous section we evaluated the polarization observed in a direction n, Ap(n, k), induced 
by just one single Fourier mode of the gravitational fluctuations. We were able to choose 
Stokes parameters Q = Ap and U — by the choice of the z axis along the direction of k, 
and defining the Stokes parameters through projections of the CMB intensity along 9 and 
0. But we can not do the same for all wavevectors. Now we must fix, for each direction 
of observation, a pair of orthogonal axis, the same for all wavevectors, to project the CMB 
intensity and evaluate the Stokes parameters. After rotation of the fc-dependent basis to the 
fixed one, for each direction of observation, there is a relatively simple relation between Q 



and U with Ap as calculated in the previous section ||11|| . Let h be a direction of observation 
on the sky. The Stokes parameters induced by a given Fourier mode of wavevector k read, 
in terms of Ap as calculated in the previous section, as 

Q(h) = A P (n,A : )cos(20g) ; U(n) = A P (n, k) sin(2^) , (3.1) 

where 0£ is, for a given direction of observation, the angle of rotation between the two basis. 

Given two directions of observation 774 and hi the polarization correlation function is 
defined as 
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C P (n h h 2 ) = (Q(nx)Q(n 2 ) + I7(ni)tf (n 2 )) (3.2) 

where (. . .) is an ensemble average. We assume that the density perturbations are gaussian. 
Thus 

(A* P (n u k)A P (n 2 , k')) = |/5(A;)| 2 (1 - - nl)e ikr ^-^8 3 (k - k') (3.3) 

where /ii and fi 2 are the cosines of the angles formed by k with -hi and n 2 respectively, 
r = r — r D is the distance to the last scattering surface, and we have written the result in 
terms of the quantity (3{k) as defined in eq. (|2.16| ). 
The correlation function can be written as 

C P (n u n 2 ) = J c?A;|/3(*)| a (l-/i?)(l-/^ e ^^)oos[2(^-^] (3.4) 
The multipole coefficients of the polarization correlation function are given by 

( a ?) = £ / / ^i^2^(ni)y Zm (n 2 )C P (ni,n 2 ) (3.5) 

In the previous section we have seen that the most interesting structure in the polarization 
of the CMB occurs for wavelengths comparable or smaller than the horizon at decoupling. 
These wavelengths subtend an angle in the sky of less than a few degrees. Small wavelengths 
have an effect on high multipoles of the correlation function expansion, / ^> 1. Thus, we 
will make simplifying approximations, valid for high multipoles only, the most interesting at 
any rate. Recall that and 2 £ are the angles of rotation from the basis used to define the 
Stokes parameters when we took the z axis in the direction of k, to the fixed basis we use 
now. They are different for each direction of observation. However, if hi and h 2 form a small 
angle 9, then cos[2(0 1 ^ — 4> 2 %j\ ~ 1 — 0(6> 4 ). Thus, replacing cos[2(0 1 ^ — <fi 2 %)] by unity in 
equation (|3.4j ) constitutes a very good approximation for high /. Under this approximation 
the multipoles have a simple expression 

(a?) « j{2l + 1) J d 3 k\A P (k)\ 2 [c l+2jl+2 (kr) + 2 ClJl (kr) + c^ 2jl „ 2 (kr)} 2 (3.6) 

where 

4(Z + !)(/ + 2) A{f + 1-1) Ajl - 1)1 

01+2 (2/ + l)(2/ + 3) ' Q (2Z-l)(2Z + 3) ' Q " 2 (2Z + 1)(2/-1) ' 1 ' 

Notice that for large /, q +2 , q_ 2 and q all tend to unity. 

Another important quantity is the power spectrum of the total polarization, W(k) defined 

as 

roc rlh 

C P (9 = 0) = (Q 2 + U 2 )= —W{k) (3.8) 

Jo k 

W(k) measures the contribution of each wavevector k to the total squared polarization. 
From equation (|3.4j) we see that it is given by 



W(k) =27rk 3 \f3(k)\ 2 J^d^l-fi 2 ) 2 (3.9) 
The total degree of linear polarization is P — (J dkWik)jk) x l 2 . 
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IV. POLARIZATION IN A COBE-NORMALIZED ft = 1 CDM MODEL 

In this section we specialize the evaluation of the polarization power spectrum and its 
correlation function multipoles to a specific cosmological model, namely a spatially-flat, Qq = 
1 cold dark matter cosmology with scale invariant, adiabatic scalar fluctuations normalized 
to the large angle anisotropy measured by COBE, under the assumption that the baryonic 
energy density is much smaller than the radiation energy density (R <C 1) at the time of 
decoupling. We will also compare our analytic results with previous numerical calculations. 

The main ingredient in eq. (|2.16|) necessary to evaluate Ap(n, k) is the value of the 
dipole in the temperature anisotropy around decoupling, which in turn can be derived from 
the monopole through the relation 

ikA T1 = -(A T0 + $) . (4.1) 

The monopole itself satisfies the equation 

x -R : k 2 . » R ■ k 2 ^ 

A T o H A T0 + — -A ro = -$ $ • (4.2) 

1 + R 3{1 + R) 1 + R 3 v ; 

In models with low baryon content, such that < 1 at the time of recombination, a 
WKB approximation to the exact solution of eq. (|4.2|) is a good approximation for most 
wavelengths of interest [2Cfl. If we define 



A = A T0 + $ (4.3) 
we get for A , in the WKB approximation, 

A = acos(co> r) + 6sin(co> r) + 2cu [ dr' sin[co> (r - r / )]$(r / ) , (4.4) 

J o 

with ujq = kc s ~ k/\f?>. We have neglected anisotropic stresses, so that $ = The 
constants a and b are fixed by the initial conditions. For adiabatic fluctuations a = |$(0) 
and 6 = 0. In order to calculate the integral above we need to know the time dependence 
of the gravitational potential. Rather than resorting to a full numerical approach, we will 
make some simplifying analytic approximations. During matter domination the potential $ 
remains constant on all scales, and the integral in (|4.4|) is easily performed. For wavelengths 



which were outside the horizon at the time of matter-radiation equality the behaviour of 
the potential is also very simple: it remains constant during both radiation and matter 
domination, but the value during matter domination is 9/10 times the value during during 
radiation domination. In this case: 

A (r) = -^$(0)cos(^ r) + ^$(0) if k -C k eq (4.5) 

For wavelengths which entered the horizon before matter-radiation equality, the situation is 
different, since the potential $ oscillates and decays. For times well before matter-radiation 
equilibrium there is an analytic solution for the potential: |30| 



$ = 3$(0) Sm(a;or) 7^ COS ^ r) (4.6) 
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If k ^> k eq the contribution of the integral is negligible, and then 

A Q (r) = -^j$(0) cos(^ot) if k^>k eg (4.7) 

With the solutions above for A (t) in the two different regimes, we can evaluate the tem- 
perature dipole from jfcAn = Ao and then use (|2.16 ) to find the polarization produced in 
this specific model by one single Fourier mode: 

f ^$(0)(0.51cJoAr D )sin(cj o r D )e ifc ^ TO - rD ) if k < k eq 
A P (A;)= ' (4.8) 

[ i§$(0)(0.51cj Ar D )sin(cJor D )e ifc ^ T °- rD ) if k > k eq 

For intermediate wavelengths the solution is not simple, so we find an approximate solution 
using eq. (|4.6|) during the radiation dominated epoch, and matching the potential to a 
constant during matter domination. 

We also have to take into account the damping of the temperature fluctuations and 
polarization at small wavelengths due to photon diffusion (Silk damping) and due to the finite 
width of the last scattering surface. We have already calculated the effect of Silk damping 
at times around decoupling. The exponential damping factor, 7 = (fc/fco) 2 is given by eq. 
( P-19Q . Since k D is time dependent, and varies very rapidly during recombination because 
of the fast change in the ionization fraction, its effect on the damping of the polarization 
should be averaged over the width of the last scattering surface. The average net damping 
factor (7) should be approximately given by 

e -<7> = / dr^e- K e { -^ ] (4.9) 
Jo dr 

The other source of damping comes from the finite width of the last scattering surface. The 
oscillations of the imaginary exponential in equation ( [2.10|) produce a cancellation in the 
integral for wavelengths smaller than Atd- To take this effect into account, the phase of 
the perturbation should be averaged over the width of the last scattering surface, and so we 
should replace e tk ^ T by 

T °dr$Vv^ r (4.10) 
dr 

— * 

This factor depends on //, the cosine of the angle between the wavevector k and the direction 
of observation n. Photons moving in a direction perpendicular to k do not suffer the damp- 
ing due to the finite thickness of the last scattering surface. To simplify the calculations, 
we average the damping factor over the angles, before performing the multipole expansion. 
Our results for small wavelengths will be qualitative, and only approximately correct quan- 
titatively. We denote by /(/cAtd) the averaged damping factor due to the finite width of 
the last scattering surface. It is given by 

/l /"To rj ft 

dfi / dr— e -V^(T-T ) (4.11) 
-1 Jo dr 

To summarize, our final expression for the polarization produced by one single Fourier mode 
is given by the product of the factor f(kAr) times expression ( |4.8|) . The damping factors 
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and the behaviour at intermediate wavelengths do not have simple analytic expressions, so 
we handle them numerically. 

Using the result for Ap, we evaluate the polarization power spectrum W(k) as given by 
eq. ( |3.9|) . Figure 1 shows the result, with the gravitational potential normalized to adjust 
the COBE-DMR measurement of the quadrupole in the temperature anisotropy, under the 
assumption of scale invariance for the scalar fluctuations. To normalize we have taken on 
large scales ^$(0) = $(r D ) = Ak~ 3 with (\A) 2 = 4!(4 2 )/5(4vr) 2 , and (a 2 T2 ) = 2 x 1(T 5 
the quadrupole in the temperature anisotropy measured by COBE. Figure 1 displays the 
main features of the CMB polarization. Very long wavelengths contribute little to the final 
polarization because of the relatively short width of the last scattering surface and the very 
tight coupling between photons and electrons prior to decoupling. For intermediate scales 
there are oscillations which follow the oscillations of the temperature anisotropy dipole, 
since those are the ones that by "free streaming" during the decoupling transition originate 
the quadrupole anisotropy that gives rise to polarization. Finally, the contribution to the 
polarization by smaller scales decays due to Silk damping and due to the cancellations 
produced by the finite width of the last scattering surface. 

We have also calculated the polarization correlation function multipoles according to eq. 
( |3.6D . Since the largest contribution and the most interesting structure in the polarization 
power spectrum occurs at intermediate wavelengths, which correspond to / 3> 1, we can 
approximate q+2, q and q_2 by unity, their large I limit. In that case the combination of 
Bessel functions appearing in eq. ( |3.6| ) can be approximated by {^:) 2 ji(kr). On the other 



hand due to the large value of r, this Bessel functions wildly oscillate, and it is a good 
approximation to replace them by their approximate average 

,. 2/ u f \2x(x 2 -I 2 ) 1 ' 2 }- 1 if x>l , AinS 
0f»> - { \f x Jt (4-12) 

The result for the multipoles of the polarization correlation function is plotted in Figure 2. 
The main features of this plot are easily understood from the shape of the power spectrum 
in Figure 1. 

The results in Figures 1 and 2 are comparable to the numerical results obtained by 
Bond and Efstathiou in reference llfl , more specifically to their figures (4b) and (7a). The 



qualitative structure of peaks is well reproduced. Our quantitative results for the heights 
and locations of the first peaks are relatively good, but accurate only up to 20 - 30 %. For 
instance, the first four peaks in our Fig. 1 appear at k « 24, 56, 100, 125 Gpc -1 respectively, 
with heights W(k) m 0.1,0.6,1.35,1.6 x 10~ n , while in Fig. (4b) of ref. JTTJ they appear 
at k « 16, 50, 80, 100 Gpc -1 with heights (after adjustment by a factor 15.4 to account 
for different conventions and normalization [TO) W(k) ~ 0.1,0.6,1.6,2.6 x 10~ n . The 
first four peaks in our Fig. 2 appear at I ~ 180,420,710,1000 respectively, with l 2 Ci ~ 



0.4,2.0,4.0,4.5 x 10 -11 , while in Fig. (7a) of ref. [|Xlj] they appear at I « 130,320,560,800 



with heights (after adjustment by a factor 7.7 to account for different normalization) l 2 Ci 
0.3,1.5,4.9,6.2 x 10 -11 . 

The quantitative discrepancies in the location of the peaks can be partially attributed 
to the simplifying approximation made in this section that the photon-baryon sound speed 
was constant. This quantity varies by around 15% from the big-bang to the time of de- 
coupling. Neglect of its time-dependence induces an error in our calculation of the phase 
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of the oscillations and thus in the location of the peaks. Besides, we cannot expect a good 
quantitative agreement of the peak heights on scales comparable to the horizon at the time 
of matter-radiation equality due to our crude approximation for the potential at these wave- 
lenghths. We can not expect good quantitative agreement for small wavelengths either 
given the rough approximations we made in this section to average the damping effects. 
Our analytic approach is not appropriate to obtain accurate quantitative predictions: due 
to our simplifying approximations we can not expect better than 20-30 % accuracy in our 
estimates of the power spectrum. The most interesting structure of peaks, however, is well 
reproduced by our analytic results, and shows that our approach can be reliably used to 
analyse the physical mechanisms that lead to polarization of the CMB, and their dependence 
upon cosmological parameters. 



V. POLARIZATION OF THE CMB IN OPEN UNIVERSES 

We have seen in Section II that the polarization of the CMB is produced during the 
process of decoupling of matter and radiation, and is proportional to the width of the last 
scattering surface Atd, and to the value of the dipole in the temperature anisotropy at the 
time of decoupling, which in turn is determined by the value of the gravitational potential 
$(td). The anisotropy in the CMB temperature also depends on these quantities, but 
differently. For instance, it is very insensitive on large angular scales to the values of td and 
At£>. It is the aim of this section to show that because of this different dependence, the ratio 
between anisotropy and polarization at relatively low multipoles is very sensitive, in an open 
universe, to the value of the matter density Qq. Besides, while the polarization of the CMB 
is produced during the decoupling transition, the present temperature anisotropy receives 
contributions not only from the Sachs- Wolfe effect, which is proportional to Q(td), but also 
from the integrated Sachs Wolfe effect (ISW), due to the time dependence of the gravitational 
potentials. In an open universe, or in a flat universe with a cosmological constant A, the 
potentials depend on time, and so the ISW effect gives an additional contribution to the 
anisotropy. We will show that this also makes the ratio between anisotropy and polarization 
dependent, upon Q . 

We shall evaluate the ratio between the multipoles of the temperature fluctuations and 
of the polarization correlation functions, and investigate its dependence upon Q for Q < 1 
in two special cases: when there is no cosmological constant, and when Qq + Q\ = 1, with 
Q\ = A/3Hq. In both cases we shall assume a scale invariant spectrum of adiabatic density 
fluctuations. 



'cum 



It has been shown [[IT|] that for multipoles such that l curv < I < Id, where l c 

irJ (1 — fi )/^o and Id = t/td, with r the distance to the last scattering surface, the 
coefficients of the multipole expansion of the temperature fluctuation correlation function 
can be approximated by 

C? = A(l + 9 -^)ir (5.1) 

where A is the normalization of the scale-invariant energy-density power spectrum of fluc- 
tuations, 3>(t£>) = Ak~ 3 , and If = [9tcI(1 + l)]" 1 , both independent of Q . The dependence 
upon Q through the function g(flo) originates in the ISW, and is given by 
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'dF\ 



2 



g(Q ) = 36vr ^ ^ -J (r - r)rfr (5.2) 

where -F(t) gives the time dependence of the gravitational potential, $(r) = $(r )-^j. 
The evolution function F is given by F(r) = D{r)/a [|^] with 

D=H Sw!rk> ^ 

The Hubble constant H satisfies 

H 2 = ( ao )4 a + a e(? Q ^ _ + A ^ ^ 

O Oq -|- G e g (3 3 

with K = — 1 in an open universe and K = in a spatially flat model. 

Our result of previous sections for the coefficients of the multipole expansion of the 
polarization correlation function is 



= - y rf d A; [c l+2 j l+2 (kr) + 2c Ul (kr) + Cl . 2jl . 2 (kr)Y (5.5) 

This result, which was derived in a spatially-flat universe, is also valid for multipoles I > 
l CU rv Indeed, for the wavelengths that contribute to these multipoles, the radial functions 
appropriate for an open universe at the time of decoupling are well approximated by the 
spherical Bessel functions of the spatially flat case @,^]. Besides, the evolution of the 
gravitational potentials until the time of decoupling is basically the same as in a flat universe 
provided that £lo z D » 1. We shall thus work under this assumption, and consider mutipoles 
such that I > l C urv only. We shall also restrict our attention to multipoles I « Id, which 
correspond to those before the first peak in Figure 2, just to be able to work analytically. 
The wavevectors that significantly contribute to these multipoles are such that kr^ << 1, 
and for them we can approximate 

\(3(k)\ 2 = [0.17$(r D HAr D sm(u T D )] 2 « [0.17$(r D V 2 Ar D r D ] 2 (5.6) 

and the multipoles can be rewritten as 

Cr = A( , ATDT ° ) 2 I? (5.7) 

' \(l + R)r 2 J 1 K ' 

with 

If = J dxx[ci +2 ji +2 (x) + 2ciji(x) + c^ 2 ji- 2 (x)] 2 (5.8) 

Thus, the ratio between the temperature fluctuation and the polarization correlation func- 
tions multipoles can be written as 



15 



with Bi independent of Qq and H . 

We want to find the explicit dependence of this ratio upon Q . The distance to the last 
scattering surface, r, is given by 

° ^ (5.10) 

TTT^ aCi 

When there is no cosmological constant this reduces to 



= 2n z D + (2^o - A)y/Qoz D + 1-1 
r H a nl(l + zn) [ ' J 

When VLqZd » 1 it can be approximated by r « 2/Q aoH . 

The dependence of and Ar^ upon f2 can De easily found exploiting the fact that z D 
and Azd are approximately independent of f2 and Taking this into account, for 

VLqZd » 1 we get 

^ D a ^o(l + ^) 3 / 2 ' a ^ /2 i/o(l + ^) 1 / 2 

Finally, the function <?(fi), which can be shown to be function of the combination (1 — 
Qq)/Qq only, is well approximated by 

^flo) « 4.87^^) (5.13) 
Using the above expressions we get for the ratio of multipoles 



Cf _nr^fn n-mi , 4.87(1 -0 )M 10 - 2 



^ = iV^(O ,0 = jV,[l + i W (5-14) 



where iV; is a normalization factor, defined so that G(Qq, I) measures the ratio between the 
temperature and polarization correlation function multipoles normalized to the value of the 
same ratio when Qq = 1. 

We plot G(flo, I) in Figure 3 for I = 30. It is clear that the dependence upon fl is very 
significant. The ratio changes by a factor of order 40 for flo ~ 0.2 and 250 for Qq ~ 0.1. 

We now repeat the calculation for a model such that Qq + = 1. g(Qo) can be 
approximated in this case by g(fto) = 0.33[(1 — f2o)/^o] 2 ' 23 - There is no simple analytic 
approximation for r in this case, but in can be evaluated numerically and, when QqZd » 1 
the result is well fitted by r m 2/Hoao^o' 39 ■ Finally the ratio of anisotropy to polarization 
multipoles is given in this case by 

S = NlG{ n , i) = Nl [i + Q-33[(i-^o)M] 2 - 23 ^o, 2 (515) 

which has a much weaker dependence upon Q than in the open universe. The result is 
plotted in Figure 4 for a fixed value of f2& and I = 30. 



16 



VI. CONCLUSIONS 



We have performed an approximate analytic evaluation of the polarization induced in 
the CMB on a wide range of angular scales by Thomson scattering prior to decoupling in the 
presence of density perturbations with adiabatic initial conditions, in a model with standard 
recombination. Eq. (|2.16|) gives the polarization induced by one single Fourier mode of the 
density perturbations, down to scales such that kAr^, rs 1. On smaller scales the finite 
width of the last scattering surface and photon diffusion damp the polarization of the CMB. 
The exponential damping factor due to photon diffusion (Silk damping), the same for the 
anisotropy and the polarization, is given by eq. ( |2.19p . We stress the fact, already pointed 
out in Ref . |[{| , that the damping factor in the anisotropy that one would derive neglecting 
the polarization dependence of Thomson scattering would be slightly incorrect: instead of 
the factor 16/15 in eq. ( [2.19| ) one would get a factor 4/5 [j37p0| . In conclusion, accurate 



calculations of the anisotropy on scales where the width of the last scattering surface is 
relevant should always include the polarization dependence of Thomson scattering. 

Eq. ( p.!6| ) displays the fact that today's polarization of the CMB was basically produced 
at the time when the tight coupling which kept the photon distribution isotropic in the 
electrons' rest frame started to break down during recombination. The degree of polarization 
is proportional to the quadrupole in the temperature anisotropy around decoupling, which in 
turn is mainly due to the "free streaming" of the dipole during the last few scatterings. This 
is manifested in eq. ( |2.16| ) through the proportionality of the polarization upon the dipole in 
the temperature anisotropy at decoupling and upon the width of the last scattering surface. 
The faster decoupling occurs, smaller the degree of polarization induced in the CMB. 

Figures 1 and 2 display our results for the CMB polarization power spectrum and cor- 
relation function multipoles respectively, for a Qq = 1 CDM model with adiabatic, scale- 
invariant scalar fluctuations normalized to the COBE-DMR measurement of the large angle 
CMB temperature anisotropy. We have assumed a standard recombination history and that 
the baryon density was much smaller than the radiation density around decoupling. The 
agreement in the height and position of the first peaks with respect to previous numerical 
computations |TT[ is very good. 

If the ionization history was as assumed here, with a tight coupling between CMB photons 
and electrons until last scattering, and no later reionization, then the relative locations of the 
peaks in the temperature anisotropy and polarization correlation functions are in a simple 
and definite relation, as we discussed in section (II. B). Thus, measurement of the relative 
locations of the peaks in the anisotropy and polarization may give additional clues to the 
ionization history. Let us also stress again that the relative heights of the peaks in the 
polarization are less dependent on parameters such as Qf, than those of the temperature 
fluctuations. Indeed, depending on the value of fib it is possible that the heights of even- 
numbered peaks of the anisotropy are suppressed with respect to odd ones, due to a partial 
cancellation between the adiabatic oscillations and the Sachs- Wolfe effect. This does not 
happen to the degree of polarization, since it is not affected by the photons' redshift. The 
absolute height of the peaks in the polarization, on the other hand, is weakly dependent 
upon Qf,; they decrease with increasing baryonic density. 

At last, but not at least, we have shown that the ratio between anisotropy and polar- 
ization multipoles is very sensitive to the value of Q in an open Universe, as evidenced in 
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Figure 3. This strong dependence upon Qq arises because anisotropy and polarization have 
a different scale- dependence. Roughly, the polarization induced by a Fourier mode scales 
as k 2 Ar D r D with respect to the anisotropy induced by the same mode, on relatively large 
scales. In an open Universe with a scale invariant spectrum of density fluctuations, and just 
for geometric reasons, a change in the spatial curvature significantly changes the relation 
between the wavenumber k of the density fluctuations and the value of the multipole / to 
which that wavenumber contributes the most, significantly changing the value of the ratio 
between anisotropy and polarization multipoles. The ratio of anisotropy and polarization 
multipoles also dependeds upon Q in an open Universe due to the time dependence of the 
gravitational potentials, which affects anisotropy through the integrated Sachs- Wolfe effect, 
while it does not affect the polarization. 
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FIGURE CAPTIONS 

Figure 1: Polarization power spectrum, as defined by eq. ( |3.9|) , normalized to the COBE- 
DMR measurement of the quadrupole temperature anisotropy, for a cold dark matter model 
with n = l,h = 0.75 and tl b = 0.03. 

Figure 2: Polarization correlation function multipoles normalized to COBE for a cold dark 
matter model with Qq = 1, h = 0.75 and f2& = 0.03. 

Figure 3: Ratio of the anisotropy to the polarization multipoles for / = 30 as a function of 
Q in an open universe with no cosmological constant, normalized to its value when Q = 1. 
Figure 4: Ratio of the anisotropy to the polarization multipoles for I = 30 multipoles as a 
function of Qq, normalized its value when Qq = 1, in a spatially-flat model with fio + ^A = !• 
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